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Thermodynamic Extrapolation of 

Rocket Performance Parameters FRANK J. ZELEZNIK 2 

NASA Lewis Research Center 
Cleveland, Ohio 


First and second partial derivatives of some thermodynamically defined rocket performance 
parameters are presented- These derivatives are useful for the extrapolation and interpolation of 
thermodynamic performance calculations for changes in combustion pressure, combustion en- 
thalpy, pressure ratio, and area ratio. The accuracy of the extrapolation or interpolation is indi- 
cated for a typical propellant combination. 


T70R SOME propellant systems, performance data are 
Jl needed over a wide range of conditions such as chamber 
pressure, pressure ratio or area ratio, oxidant to fuel ratio, 
and, occasionally, initial enthalpy. The thermodynamic cal- 
culation of rocket performance is usually sufficiently difficult 
to require the use of automatic computers. Because of the 
cost and time involved in the computations, it is often not 
feasible to explore the complete range of independent vari- 
ables. For this reason, the problem of extrapolating or inter- 
polating a limited number of performance calculations occurs 
frequently. 

The use of partial derivatives in extrapolating performance 
data with considerable accuracy for moderate changes in the 
independent variables was discussed previously by Gordon 
and Huff (l). 3 In Ref. 1, expressions were derived for first 
partial derivatives of the logarithm of specific impulse, charac- 
teristic velocity, and area ratio with respect to the logarithm 
of chamber pressure. In addition, first- and second-order 
corrections to impulse for a change in combustion enthalpy 
were also given. A similar first-order correction for the effect 
of combustion enthalpy (or heat of formation) on impulse was 
described by Gordon (2) . 

In this report, the list of partial derivatives is extended to 
include the derivatives of specific impulse 7, specific impulse 
in vacuum /vac, area per unit mass flow rate A ( w , area ratio e, 
and characteristic velocity c* with respect to combustion 
pressure P c , combustion enthalpy h C} pressure ratio P c /P, and 
area ratio e. Examples are given to illustrate the use of the 
derivatives. 
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Analytical expressions for the derivatives of parameters 
with respect to oxidant to fuel ratio 0/F can be derived; how- 
ever, the method of derivation of these expressions is some- 
what different from that used to obtain derivatives for a fixed 
O/F, and, therefore wall not be presented in this paper. All 
derivatives in this paper are for a fixed O/F. 

Effect of Change in Initial Conditions 
on End Point of a Process 

All the thermodynamic properties of a system of known 
composition can be specified uniquely in terms of any two 
thermodynamic functions, say a and 0. Thus the system can 
be represented by a point in a two-dimensional space with 
coordinates (a,jS). At any point in (a,fi) space, not only are 
all the thermodynamic properties of the system known, but 
it is also possible to determine the rate of change of these 
properties along some curve in (a, ft) space. Thus if ^ is a 
third thermodynamic function, the derivative (d ft /da) 4 , ex- 
presses the rate of change of with respect to a change in a 
along a curve of constant \f/. This partial derivative is the 
usual thermodynamic first partial derivative discussed in 
thermodynamic textbooks. 

Thermodynamic Processes 

A process in thermodynamics means that a system orig- 
inally at some point (ao,/3o) has been moved to a new point 
(a,/3). An infinitesimal process can be completely charac- 
terized by a derivative of the form (d/3/da)^. A finite process 
can be specified by giving a starting point (ao,0o), a P a *h> 
say a curve of constant and one of the coordinates of the 
end point a. For a given path, the only independent vari- 
ables of the process are the coordinates of the initial point 
(a Ot 0o) and a coordinate, say a, of the final point. The other 
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coordinate of the final point, ft and all other thermodynamic 
functions are dependent variables in the process. 

Let be any dependent thermodynamic variable associated 
with the end point of a finite process. The rates of change of 
<p with respect to the independent thermodynamic variables 
of the process are of two different types. For a change in a , 
the usual type of derivative (b(p/da)$ is needed. However, 
for a change in one of the coordinates of the initial points, 
derivatives of the type (d<p/dao)p 0 and (d(p/dj3 0 ) a(j are needed. 
Expressions for the latter type of derivative are obtained in 
the following section. 


Derivatives Associated With a Finite Process 


As indicated in the previous section, the process is defined 
as taking place along a curve of constant \p. Therefore 

Mao, ft) = iKa, P) = <p) [1] 

Whatever change in results from a change in the initial 
point (a 0 ,ft) must be equal to the change in \p. This may be 
expressed in differential form as 

(Is). + (Is ). m - (HO, 4 “ + (H)."" 121 

Imposing in turn the conditions of constant a 0 and constant ft 
gives the following desired expressions for the partial deriva- 
tives of a function at the end point of a process with respect to 
the initial coordinates : 


/<V\ = W/Wa (d<p\ (da\ 

\M/« o ^ W* \^ft/ ao L J 

\ = Wo/d<Xo)3o (d<p\ (da \ 

\daoJf} o (d^/dip) a Vdao/jSo 

In Eqs. [3 and 4] all the derivatives except (da/dft) ao and 
(da/da 0 )/ 3 0 are the standard thermodynamic first derivatives 
and can be immediately evaluated. The two exceptions can 
be evaluated for a specified form of the relation 


a = a(a 0 , ft) [5] 

Two forms of Eq. [5] are considered in this paper: a = k x 
and a — k 2 a^ where hi and k 2 are constants. For a = k h Eqs. 
[3 and 4] reduce to 



For the particular choice a = k 2 a^ Eqs. [3 and 4] give 



It may be seen that the right-hand sides of Eqs. [6 and 8] 
are identical. 


Application to Rocket Performance 

First Derivatives of Thermodynamic Functions With 
Respect to P c , h c , and PJP for an Isentropic Process 

In this paper, Eqs. [6-9] are applied to an isentropic proc- 
ess. For an isentropic expansion from P c and h c to an exit 
pressure P, the following correspondences to functions in 
Eqs. [6-9] apply: s c = \[/ 0 , P c = ao, P = a, h c = ft, and 
s = \p. From Eqs. [6 and 8] 


\bh c ) p c ,p \dh c J p e ,p c /p 

_ /5sc\ //bs\ = 1 

\bhjpj \d<p) P T c (ds/d<p) P 

FromEq. [7] 

/d<p\ = (ds e \ //ds\ _ Zh 

Wc) he ,p WJ hc / W/p / Wa 

From Eq. [9] 



[ 10 ] 


111] 


[ 12 ] 


For any <p selected, the right-hand sides of Eqs. [10-12] 
are in terms of the usual first partial derivatives. Bridgman 
(3) presents a convenient scheme (given in many thermody- 
namic textbooks) for expressing all first partial derivatives in 
terms of three first partial derivatives, namely, ( dh/dT)p = 
c p , (dv/dT) P , and (dv/dP) T - An equation of state is required 
to evaluate these derivatives. In rocket performance calcu- 
lations it is convenient to assume that the reaction products 
obey an ideal equation of state with a variable molecular 
weight 

Pv = RT/M [13] 


It has been found useful in performance calculations to 
calculate molecular weight derivatives rather than volume 
derivatives to indicate changes in composition. From Eq. 
[13] there follows 


(dv \ v r (d In M\ T 
\dT/p “ T L 1 ' \d In T )p\ 

(dv \ _ — v r / d In M \ T 

\dp)r “ L 1 + \a In P ) tj 


[14] 

[15] 


A further discussion of the molecular weight derivatives and a 
method for their numerical evaluation is given by Gordon et 
al. (4). With the aid of Bridgman’s tables (3) and Eqs. 
[13-15], all other first partial derivatives can be expressed 
in terms of c p (a In M/d In T)p and (a In M/d In P) T . As 
an example, assuming <p = T, Eq. [12] becomes 



In a manner similar to the previous example, derivatives 
were obtained from Eqs. [10-12] for <p = T, h, and p. Ex- 
pressions for the logarithmic form of these derivatives are 
given in Table 1. The logarithmic form was selected to pre- 
sent the results in dimensionless form. Expressions for 
molecular weight derivatives are also given in Table 1. These 
are obtained by considering molecular weight to be a function 
of temperature and pressure and expanding by the chain rule 
in differentiation. Derivatives of the form [d In <p/d In 
(Pc/P) ]PcM = In #>/(& In (Pc/P) ]. = - (d In <p/d In P) 6 are 
also given in Table 1. Henceforth, derivatives at constant 
P e and h c will be shown as derivatives at constant s. The use 
of these derivatives will be discussed in a later section. 

A word of caution is required about the derivatives with 
respect to In h e . If a reference level for enthalpy is chosen so 
that he is negative, the logarithmic derivative must be re- 
placed by a derivative with respect to h c . 
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Table 1 First derivatives of some thermodynamic properties and rocket performance parameters 
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Table 2 Second partial derivatives of 7 2 


l~ d(d/ 2 /d In P c )hc,Pc/p l T 2 7j7^\ 7 d In T \ / d In 7 \T^ 7^ ^ M< \ 7 d In M \ iT 0 

L d In Pc J kc , Pc/P ( \d In Pc ) h c , Pc/P \& In P c ) h e , Pc/P \d In P c / P )• LM C \d In P c ) h e \d In P c ) h *, p c / p \) 

1 = r ^iydh c ) PetPc/F i = 1 * t r/ d in ta _ / a in t \ j 

J Pc, Pc/P L a In P c J hc,Pc/P h c - h T c L \a In P c )hc \d In P c )hc,Pe/PJ 

1 = r a(a/yainP c /P). n _ gf / a in / \ r/ a in r \ __ / a in jf \ J 

J, L b \ nP c J hc . Pc/p \a In ( P c / P)/ M L \d In P c ) hc , Pc/P \<> In P c ) k c ,p c /p] 


[" ^ d/2 / d * n P c)hc ,P c /P~ 

L a^c 

r a(a/ya in P c )he,p e / p~ 

L a In P,/P 

r a(a/ya^)F Ct p c /p ~[ = /» _r / jl^ _ i\ 

L a/ 2 c JPe^c/P X c - h T c 2 \C/> e C P ) 


7 2 P7 7 


r a(a/ 2 /a/t c )p c> p c /p ~| = r a(a/ 2 /a in p c /p), ~| 

L a In P c /P J, L afc e Jp^.iVP h c ~ h Mc t 

r a(a/ya in p c /p), 1 _ 2J / i - 7 \/ a in/ \ 

L a In p c /p J, V 7 / \a In p c /p ) , 


r r /a in a/ \ * 
p c L \a in r ) p _ 


« (b In M c \ 

2R fi 

( b In Me 

\dln P c )hc 

M c cp c L 

\ d In Tc t 

6 /din Tc\ 


/d In Mc\ 

\d In Pc) he 

- — i 

M'Cp. 

\bln Tc ) P e 


7 c - 1 


M c c Vq 


Derivatives of Rocket Performance Parameters With 
Respect to P c , h c , and P c /P for Isen tropic Expansion 

Parameters that are usually of interest in performance cal- 
culations are defined in the following equations, where it is 
assumed that a consistent set of units is employed : 


/ = [2 (h e -h)]'t* [17] 

(A/w) = 1/pu = 1 /pi [18] 

e = {A/w)/(A/w), = A/ A t [19] 

c* = PUMt [20] 


^’C-SOC-S- 7 )] ™ 

Detailed expressions for the second partial derivatives of 7 2 
obtained by means of [24b] are given in Table 2. 

The expressions for the first derivatives of the remainder of 
the performance parameters in Table 1, unlike specific im- 
pulse, include first derivatives of thermodynamic functions. 
The second derivatives of these performance parameters are 
therefore not included in this paper, inasmuch as they involve 
second derivatives of thermodynamic functions which are 
generally not available. 


/vac = / + (A/W)P 

C F = I/c* 


[ 21 ] 

[ 22 ] 


First Derivatives of Rocket Performance Parameters With 
Respect to P c , h c> and e 


A specific set of units and associated dimensional constants for 
Eqs. [17-22] is given in Ref. 4. 

Derivatives of the logarithm of 7 can be obtained directly 
from the derivatives of the logarithm of h in Table 1 . Thus, 
for example 


/ d In 7 \ = —h / bin h \ 

\d In Pc)p c /p,h c 7 2 \d In P c )p c /p,h c 


[23] 


Eqs. [18-22] can be used to obtain derivatives of the re- 
maining parameters in terms of the derivatives of the loga- 
rithm 7 and the thermodynamic functions. The results are 
included in Table 1 . 


Second Derivatives of Specific Impulse With Respect to 
P c , h c , and PJP 


The previous discussion of thermodynamic processes 
showed that a process could be specified by giving the value 
of two thermodynamic functions at the initial point, a path 
connecting the final point and initial point, and one of the 
coordinates of the final point. Thus for a given process, any 
parameter of the process X can be considered a function of the 
coordinates of the initial point and the coordinate of the 
final point. For the isentropic expansion of gases in a rocket 
nozzle, it is permissible to write 

X = X(P C , h c , PJP) [25] 

and for area ratio in particular 

e - c(P„ K PJP) [26] 

Eliminating PJP between the previous two equations gives 


An examination of Table 1 shows that the expressions for 
the derivatives of the logarithm of 7 include only the thermo- 
dynamic functions h, T , and M. Therefore, the second 
derivative of the logarithm of 7 can be obtained by the 
methods previously discussed. However, for purposes of 
extrapolation, 7 2 is a better form than In 7. (This will be dis- 
cussed further in a later section.) The first and second 
derivatives of 7 2 can be expressed in terms of logarithmic 
derivatives as follows : 


^ = 2/^ 
by by 


[24a] 


X = X(P C , h c , e) [27] 

where the functional forms of Eqs. [25 and 27] are not the 
same. 

If the functional relation in Eq. [25] were knowm, then this 
would provide an alternate method of obtaining the first, 
second, and fourth columns of derivatives in Table 1, i.e., 
(d In \/b In P c )p c /p,k c . , (b In X/d In h e )p c /p,p e , and [b In \/b 
In (PJP ),] . Similarly, if the functional relation in Eq. [27] 
were known, the partial derivatives ( b In X/d In Pc) e ,h e , (d In 
X/d In h c ) e ,p c , and (d In X/d In e), could be obtained. How- 
ever, expressions for these last three derivatives need not be 
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qotained explicitly inasmuch as they can be related to the first 
set of three derivatives as shown in the next paragraph. 

The total differentials of Eqs. [25-27 ] are 


d In X = 


(^v~w) dlnP ° + (MHr) x 

\d In Pc)p c /p,k c \£> In hj p c /p,p c 

* ln *- + [dr^]/ l ”' I - / ' p |28 ' 
‘ ,ln ''- + [Mrfe)l <,l " P ' /P I29 I 


/ c) In X \ /d In X \ 

d In X = ( —j — — ) d In P e + ( — ) X 

\d In Pc/t, he \d In h c / t , Pc 

d In h c + ^ d In e [30] 

Eliminating d in e in Eq. [30] with Eq. [29] and then compar- 
ing the coefficients of the differentials with those of Eq. [28] 
with some rearrangement, the following relations give: 

/ 5 In X \ 5 In X 1 

Vd In Pc),, he ~ \2> In Pc)p c /p,h c ~ Id In (P./P) J. X 

(d In e/d In P c )p c /p,h c , , 
[d in t/d In (P e /P)]. 1 J 
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H 2 (g) + l/30 3 (g) 
{h c = 629.1 CAL/G, 
P c =IOO PSIA)-^, 


O REFERENCE POINT 

ACCURATE / 

EXTRAPOLATED/ 


'-H 2 (g)+ 1/2 0 2 (g) 
(h c = 0 CAL/G, 
P c =IOOO PSIA) 


Table 5 Comparison of accurately calculated characteristic 
velocity data for stoichiometric H 2 -0 2 with data extrapolated from 
reference point 

h c — 0 cal/g h c = —190.6 cal /g h c = 629.1 cal/g 


Pc, 

psia 

Accu- 

rate 

Extrapo- 

lated 

Accu- 

rate 

Extrapo- 

lated 

Accu- 

rate 

Extrapo- 

lated 

1000 

7208 a 

7208 

7075 

7077 

7625 

7641 

600 

7165 

7166 

7035 

7034 

7573 

7599 

100 

7008 

7017 

6887 

6886 

7386 

7451 


a Reference point for extrapolation. 


''H 2 ( i ) + 1/2 o 2 u> 

(h c = - 190.6 CAL/G, 
P c =IOO PSIA) 


1-1 .lil.l L_j-J_l1i1i 1 I ■ I ■ hill 

too 1000 10,000 

PRESSURE RATIO. P c /P 


Fig. 1 Comparison of accurately calculated specific impulse data 
for stoichiometric H 2 -0 2 or H 2 -0 3 with data extrapolated from the 
reference point 


/dlnX\ _ / d In X X __ / 

\d In h c ) t,p c \d In hj p c /p,p e \ 


/ d In X \ 
\d In e) 8 


/ dlnX \ 
e/P , Pe ~ \i> In (Pc/P) /. X 
(d In e/d In h c )p c /p,p c 
[d In e/d In (P./P)}, 

[d In X/d ln(P e /P) ], 

[a In e/ In (P e /P)]. 


Second Derivative of Impulse With Respect to Area Ratio 

An approximate second derivative of impulse with respect 
to area ratio for constant P c and h c can be obtained by dif- 
ferentiation of [33], assuming that [5 In c/d In (P c /P)] 8 is 
constant. This is a very good approximation for pressure 
ratios greater than 10. With this assumption, there follows 


" d 2 (7 2 ) 1 

_d (In e)*l 


( 1 - y \[ dlnI 1 

V 7 Jld\n(P c /P)[ 
d In c I 2 

_d In {P c /P) J, 


Partial Derivatives for Extrapolation 

The first and second partial derivatives given in Tables 1 
and 2 can be used for extrapolation purposes by assuming 
that the change in any one of the parameters X or functions <p 
can be obtained from a truncated Taylor series. Thus, for 
any function /, neglecting derivatives higher than the second 


ffe) = /fe°) 4- 


1 1~)° (*. - 

! t W <* 

2 i,j \dXiC)Xi/ 


(Xi - xViXi - Xj°) [35] 


When only first derivatives are available, the last term in Eq. 
[35] is omitted. 

The accuracy obtained with Eq. [35] depends considerably 
on the form of / and Xi. Those forms whose first derivatives 
are the most nearly constant over the range of interest can be 
expected to give the most accurately extrapolated values. 

An indication of the desirable forms for / and Xi can be ob- 
tained by an examination of some performance data and 
derivatives. In Table 3, data are given for stoichiometric 
H2{g)-0 2 (g) for P c = 1000 psia assuming equilibrium com- 
position during isentropic expansion. Table 3 is the direct 
output from an IBM 704 program. 

The symbols used in Table 3 for some parameters are some- 
what different from those used in the rest of this paper be- 
cause the IBM printer does not contain characters such as 
lower case letters, Greek letters, subscripts, or superscripts. 
The following examples illustrate the differences: 

(DLI/DLPC)PC/P = (d\nI/d\nP c ) Pc/Pthc 

(DLCS/DHC) PC/ P - (din c*/bh c ) Pc /p tPe 

(DLAR/DLPCP)S = [d In e/d In (P c /P)] s 

It can be seen from Table 3 that the derivatives [d In e/d In 
(P e /P)] 9 are very nearly constant over a considerable range 
of pressure ratio. It is to be expected, therefore, that a good 
form of Eq. [35] for extrapolating area ratios to other pressure 
ratios is 

ln { = ° n {) ° + [5 In (/I [ ]n (f) ~ b (f)°] [36] 

In contrast, the large differences in the values for [d ln 1/ 
d ln {Pc/P) ]s indicate that using ln 7 for /and ln(P c /P) for x 
in Eq. [35] would not give particularly good results. A more 


Table 4 Comparison of accurately calculated area ratio data for stoichiometric H 2 -0 2 with data extrapolated from reference point 


Pc/P 

he 

Pc 

Accurate 

= 0.0 cal/g 
= 1000 psia 

Extrapolated 

he — 
Pc = 
Accurate 

— 190.6 cal/g 
100 psia 

Extrapolated 

he ~ — 
Pc = 
Accurate 

629.1 cal/g 
100 psia 

Extrapolated 

10.0 

2.440 

2.228 

2.468 

2.289 

2.466 

2.337 

40.83 

7.009 

6.941 

7.151 

7.130 

7.169 

7.279 

68.05 

10.50 

10.49 

10.75 

10.77 

10.81 

11.00 

100.0 

14.31° 

14.31 

14.69 

14.70 

14.81 

15.01 

300.0 

35.05 

34.75 

36.28 

35.70 

37.02 

36.45 

400.0 

44.39 

43.84 

46.04 

45.04 

47.22 

45.98 

600.0 

61.92 

60.83 

64.44 

62.49 

66.66 

63.80 

800.0 

78.38 

76.75 

81.77 

78.83 

85.21 

80.49 

1000.0 

94.08 

91.90 

98.31 

94.40 

103.2 

96.38 


0 Reference point for extrapolation. 
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nearly constant set of derivatives can be obtained in this case 
by multiplying the logarithmic derivatives by 7, i.e. 


r d/ 1 7 

r A In 7 1 

lb In (P e /P)J. 

_c> In {PJP) J. 


A still more nearly constant set of derivatives are 

[d In (PJ P)1 = 2/2 [d In (P./P)]. [38] 

For specific impulse, therefore, a good form of Eq. [ 35 ] for 
extrapolating to other pressure ratios is 

i r d 2 1 2 i° 

2fe7FTd . (Aln ^ /P)2 [39] 

A different situation occurs for temperature. In Table 3 
the derivative [d In T/b In ( P e /P)] e is fairly constant over 
part of the range of interest (from PJP = 1 to about 100), 
whereas another derivative form [bT/b In {PJP)] S is more 
nearly constant over the remainder of the pressure ratio in- 
terval. A further complication arises because the derivative 
of temperature or logarithm of temperature with respect to 
pressure ratio has a maximum. Extrapolation may, there- 
fore, give very poor results no matter which derivative form 
is used. For most chemical systems, the derivatives [d In T/b 
In {PJP) ]* are more nearly constant as a function of pressure 
ratio than are the derivatives [ bT/b In (PJP) ]*. 


Numerical Examples of Extrapolation 
Specific Impulse 

The use of the first and second derivatives permits the ex- 
trapolation of specific impulse data with considerable ac- 
curacy. This is illustrated by two numerical examples. 
Both examples start with the data in Table 3 for a pressure 
ratio of 100. 

Example 1 

Extrapolate 7 from PJP — 100 to PJP = 1000 (A In 
{PJP) = 2.30259). 

From Table 3, 7 = 382.0, [b In I/b In (PJP)] S = 0.08396, 
and 7 = 1-1215. FromEq. [24], [d7 2 /dln (PJP)] 8 = 24,504, 
whereas from the equation in Table 2, [d 2 (7 2 )/d(ln PJP) 2 ] g - 
—2654.6. From Eq. [35] 

7 2 - (382.0) 2 + 24,504 (2.30259) + i(~ 2654.6) (2.30259) 2 

7 2 - 195 309 

7 = 441.9 

This is an excellent agreement with the accurately calcu- 
lated value of 442.2 given in Table 3. 

Example 2 

Extrapolate 7 from PJP = 100 to 1000 (A In {PJP) = 
2.30259), from P c = 1000 to 100 psia (A In P c = -2.30259), 
from h c = 0 to —190.6 cal/g (from gaseous to liquid pro- 
pellant). 

Using the data in Table 3 for combustion chamber condi- 
tions and for a pressure ratio of 100, the three first derivatives 
of 7 2 from Eq. [24] are 

{bP/b In Pc) Pe /P, he - 2600.4 

{bP/bhc) p e /p,p = 28,648 

[bP/b\n(PJP)] t = 24,504 
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a Points used to obtain interpolated data. 


Numerical values for all the quantities that appear in the 
expressions for the second derivatives in Table 2 are given in 
Table 3 with the following three exceptions: 

= 0.015198 (see footnote, Table 2) 

he 

= 0.041197 (see footnote, Table 2) 

= 0.007821 (see Table 1) 

d In Pc) p c /p t h e 



The six second derivatives of 7 2 are (with subscripts omitted 
for brevity) 


d 2 (P)/d (In P c ) 2 
d 2 (7 2 )/d In PMc 
d 2 (7 2 )/d In P c b In (P c /P) 
d 2 (P)/dh c 2 
5 2 (/ 2 )/5W In (Pc/P) 
d 2 (7 2 )/d[ln (P c /P )] 2 


- 182.69 

1.5610 

162.67 

-0.0036205 

5.3690 

-2654.6 


Using the first and second derivatives in Eq. [35] gives 

P = (382.0) 2 + 2600.4(- 2.30259) + 28,648(-0.1906) + 
24,504(2.30259) + §(-182.69) (-2.30259) 2 + 
§(-0.0036205) (-190.6) 2 + §(- 2654.6) (2.30259) 2 
+ 1 .5610(— 2.30259) ( — 190.6) + 162.67 X 
(-2.30259) (2.30259) + 5.3690(- 190.6) (2.30259) 


P = 180,778 
7 = 425.2 


This is in very good agreement with an accurately calculated 
value of 426.5. 

Some additional numerical comparisons of accurate data 
with data extrapolated for various changes m P c , h c , and 
Pc/P are given in Fig. 1 for specific impulse, in Table 4 for 
area ratio, and in Table 5 for characteristic velocity. The 
extrapolated data are all obtained from one accurately cal- 
culated reference point corresponding to P c — 1000 psia, h c = 
0, and Pc/P = 100. The derivative forms used were first and 
second derivatives of P, first derivatives of In e, and first 
derivatives of c*. It is apparent that the extrapolated values 
are in excellent agreement with the accurately calculated 
values for a considerable range of extrapolation. 


Partial Derivatives for Interpolation 

Derivatives can be used to increase the accuracy of inter- 
polation in a specified range. This is because each derivative 
is approximately equivalent to having an additional point 
in the specified interval. For example, if only functions are 
known at two points, only linear interpolation is possible. 
However, if the first derivatives of these functions are also 
known at the two points, cubic interpolation is possible. 
With second derivatives also known, quintic interpolation is 
possible. Thus, for example, the following equation may be 
obtained : 

P = A + B[\n(P c /P)] + C[ln (P c /P )] 2 + 

D[ln (P e /P)Y [40] 

where the coefficients A , B, C, and D are determined by the 
solution of four simultaneous equations involving the values 
of 7 2 and first derivatives of P, each at two pressure ratios. 

To illustrate the accuracy of interpolation which can be ob- 
tained by this technique, values were interpolated for T, M, e, 
and 7 at pressure ratios of 300, 400, 600, and 800 using cubic 


equations derived from the data of Table 3 at pressure ratios 
of 100 and 1000. The functional forms for I and e which gave 
the best results as functions of the logarithm of PJP are 7 2 
and In c. Both the linear and logarithmic forms gave es- 
sentially the same results for T and M. 

The interpolated results are compared in Table 6 with the 
accurately calculated values of Table 3. Included in Table 6 
are values of C F and 7 va c Values of C F were calculated from 
c* and the interpolated values of 7 by means of Eq. [22], using 
consistent units. Values of I vac were calculated from c* 
and the interpolated values of 7 and e by means of an alternate 
form of Eq. [19]: 

7v* c = I + c*e/(P e /P) [41] 

again using consistent units. As may be seen from Table 6, 
the interpolated results are in excellent agreement with the 
accurate results. 

In the case of specific impulse, a quintic equation can be ob- 
tained from data at two pressure ratios, inasmuch as both first 
and second derivatives are available. Specific impulse 
values interpolated over the entire range of pressure ratios 
from 10 to 10,000 from a quintic formed from data at these 
two points were within 0.1 Ib-sec/lb of accurately calculated 
values. 

Data are often desired at assigned area ratios. By use of 
derivatives obtained from Eq. [33], accurately interpolated 
parameters corresponding to assigned area ratios can be ob- 
tained in a manner similar to that just discussed for assigned 
pressure ratio interpolation. 


Nomenclature 


A, B, C, D 
A/w 

— 

coefficients, Eq. [40] 

= 

area per unit mass flow rate 

c* 


characteristic velocity 

Cp 

= 

heat capacity at constant pressure per unit mass 

c F 

= 

thrust coefficient 

h 


enthalpy per unit mass 

7 


specific impulse, Eq. [17] 

7 vac 

= 

specific impulse in vacuum, Eq. [21 ] 

k lt hi 

= 

constants 

M 

= 

molecular weight 

O/F 

= 

oxidant to fuel weight ratio 

P 

= 

pressure 

Pc/P 

= 

pressure ratio 

R 

= 

universal gas constant, 1.98726 cal/mole-°K 

s 

= 

entropy per unit mass 

T 

= 

absolute temperature 

u 

= 

velocity 

V 

= 

specific volume 


- 

any variable 

a,(3,<p,rp 

= 

any thermodynamic variable 

e 

= 

area ratio 

X 

= 

process parameter 

P 

= 

density 

7 

= 

isentropic exponent (d In P/d In p) s 

Subscripts 

c 

— 

combustion chamber 

0 

= 

initial coordinate of a finite process 

t 


throat 


Superscript 

0 = reference condition 
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